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Using a discrete model of 2-dimensional easy-plane classical ferromagnets, we propose that a 
rotating magnetic field in the easy plane can switch a vortex from one polarization to the opposite 
one if the amplitude exceeds a threshold value, but the backward process does not occur. Such 
switches are indeed observed in computer simulations. 



There is a growing interest in non-equilibrium dynam- 
ics of quasi-two-dimensional magnetic materials 
Many magnetic properties of these materials are well- 
described by the classical two-dimensional Heisenberg 
model with easy-plane symmetry. In this model vortices 
play a very important role. They cause a topological 
phase transition || , and they contribute to the so-called 
central peaks in inelastic neutron scattering experiments 
that arise from the translational motion of vortices 
lOj. 

There are two types of static vortex solutions, de- 
pending on the anisotropy strength 8 (see below) of the 
Heisenberg exchange interaction in-plane vortices 

for which all spins lie in the easy xj/-plane, and out-of- 
plane vortices which exhibits a localized structure of the 
^-components of the spins around the vortex center. In 
addition to the vorticity q — ±1, ±2, . . ., the out-of-plane 
vortices have a second topological charge p. This is de- 
noted as "polarization" because its sign determines the 
side of the xy— plane to which the out-of-plane vortex 
structure points. 

The aim of the present Letter is to investigate the 
internal dynamic of the vortices in the discrete two- 
dimensional Heisenberg ferromagnet in the presence of 
an ac magnetic field in the easy-plane. It is shown that 
switching between vortex states with different polariza- 
tion occurs if the amplitude of the field is larger than a 
threshold value. 

The easy-plane Hamiltonian for classical spins Sft = 
S^sin 9ft cos sin dn sin <j>,-j, cos 9^} located on sites ft = 
{n x , n y ) of a quadratic lattice has the form 



1 - Ml Jl - M|_ s cos(<F, - $ n -a)} 



(1) 



where 8 is the anisotropy parameter (0 < 8 < 1), J > 
is the exchange constant and a is the vector which con- 



nects a site with nearest neighbors. 55 = = cos 9a is 
the on-site magnetization. In what follows we set J = 1 
and 5 = 1. 

It is known (H]|ll|] that for 8 > 8 C where the critical 
value of the anisotropy parameter 8 C depends on the lat- 
tice type(e.g. for square lattices 8 C ~ 0.3) the in-plane 
vortex with = and the azimuthal angles sat- 
isfying the equation 2s sm (^« — ^n-a) ~ 0' ^ s stable. 
For 8 < 8 C the in-plane vortex becomes unstable and 
an out-of-plane vortex is created. It exhibits a local- 
ized structure of the components around the vortex 
center. In the case of a circular system of the radius L 
with free boundary conditions the azimuthal angles <f>^ 
for both types of vortices are approximately given by 



<Pfi = q arctan — 

n r — X 



— q arctan 



n T . - X 



, (2) 



where a constant phase has been omitted. X and Y are 
the coordinates of the vortex center, and X — XL 2 /R 2 , 
Y = YL 2 /R 2 (R 2 =X 2 + Y 2 ) are the coordinates of the 
"image" -vortex. In the case of fixed boundary conditions 
the sign in front of the second term in Eq. (Q) is reversed. 

We are interested here in the vortex dynamics under 
the influence of a spatially uniform in-plane ac magnetic 
field h(t) = h(cosojt,s'mut,0). 

The interaction of the field with the spin system has 
the form 



V(t) = -hJ2J l - M H cos ( $ « - ut )- 



(3) 



The spin dynamics is described by the Landau-Lifshitz 
equation 

d t dH 



1 - M 2 d^ a 



^-^{H + V^-^l-MD—. (4) 
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The last terms in Eqs (Q) represent damping Jill]. 

To clarify the behavior of out-of-plane vortices in the 
presence of the ac field, we have numerically integrated 
the Landau-Lifshitz equation (||) for a large square lat- 
tice in which we cut out a circle with radius L = 24 
using both free and fixed boundary conditions with 5 = 
0.1, 7 = 0.002. We used relatively weak ac fields so as not 
to change the ground state significantly. The integration 
time was 12,000 time units with time step 0.01. First 
we used an out-of-plane vortex with polarization p=l as 
the initial condition and a clockwise rotating magnetic 
field with the frequency to = —0.1: This is close to the 
frequency of the lowest radially symmetric eigenmode in 
the presence of a vortex pT| , ^3[ , ^4[ . We observed that for 
all h < h cr — 0.0025 the vortex with p=l remains the sta- 
ble configuration but for h > h cr a flip to the state with 
the opposite polarization (p — — 1) occurs ( Fig 1 ). Us- 
ing the same initial condition but changing the direction 
of rotation of the magnetic field (to — 0.1) we observed 
the switching only when h > 0.02. But in contrast to the 
previous case when the vortex after switching had a well- 
defined core structure, now the out-of-plane structure is 
almost completely destroyed by spin waves. 

Another set of simulations was performed using a static 
out-of-plane vortex state with the same h but polariza- 
tion p = — 1 as initial condition. We found that the 
flip occurs only for a counter-clockwise rotating magnetic 
field to > 0. The results of extensive simulations may be 
summarized as follows for both types of boundary condi- 
tions: 

• Flips between oppositely polarized states take place 
under the action of the ac magnetic field when 
h > h cr ( Fig. 2). 

• The threshold value h cr depends on the product lo p 
and not on the vorticity. The threshold value in the 
case when the polarization vector is anti-parallel to 
the angular velocity vector Co = (0,0, co) is much 
smaller than when these vectors are parallel. 

• Flips are uni-directional. When top < 0, the final 
state is characterized by a well-defined core struc- 
ture, while for top > the core structure is de- 
stroyed. 

The basic reason for the switching can be easily under- 
stood by using the frame of reference which rotates to- 
gether with the magnetic field. In this frame there exists 
an inertial force equivalent to a magnetic field aligned 
along the angular velocity to. Then the vortex states 
with different polarization are nonequivalent and switch- 
ing processes become energetically favored. This can- 
not explain, however, why the threshold of the switching 
is a non-monotonic function of the frequency to ( Fig. 
2). To gain deeper insight we need in a reduced form of 
the Hamiltonian (fy) which takes into account both types 



of vortices: in-plane and out-of-plane. As a topological 
charge,}?, is conserved in the continuum limit only. Thus 
the switching between states with different polarization 
is due to lattice discreteness. 

We consider the near-critical case \(S— S c )/(1— S c )\ <C 1 
when the out-of-plane spin deviations are small, 
and assume also smooth dependence of the deviations 
4>ti = $n — ®ft from the static vortex structure on the 
spatial variable n. In this case, applying the transforma- 
tion 

M K = ^2Cft v m v , <fe = ^JC^Vv (5) 

where the coefficients £n,v, K-n,v satisfy the set of equa- 
tions Ch, v = 2 J2s - Kn-s,v) cos($a - $°_ s ), 
^!C H ,„ = 2E 5 (-(1 - 8)£*-it,u + £n,u cos($a - S\ 
we transform the harmonic part of the Hamiltonian (Fy) 
obtained in the vicinity of the static in-plane vortex, 

H = \ X>« - fc-zf cos($£ - - 

J2 (A M n M^ s - Mi cos(<f£ - (6) 

to the principal axis coordinates 

In Refs [ p^lJTil the linear eigenmodes of the easy-plane 
ferromagnet with Hamiltonian (Q) in the presence of a 
vortex were investigated. The lowest radially symmet- 
ric mode, which is localized near the vortex center and 
describes the in-phase motion of the core spins, becomes 
soft when S approaches S c . In other words, this mode is 
responsible for the in-plane vortex instability. 

The corresponding eigenvalue, say Hi, becomes neg- 
ative when S < 6 C : Hi = B (5 — S c ) with B a nu- 
merical coefficient. Inserting the transformation (^) into 
the Hamiltonian H± = H — Hq and keeping only the 
soft eigenmode v — 1, we obtain an effective soft-mode 
Hamiltonian 

H s = l -{i,l+^m\) + ^m\, (7) 

where terms m\il)\ and ipf which are unimportant in the 
near-critical case and all higher order terms have been 

omitted. A = | £ 3i3 (cl, - C%_^) 2 cos($« - $°_ 3 ) 

is a positive constant. 

We see from Eqs (|^) and (||) that a spatially uniform 
in-plane magnetic field cannot excite the radially sym- 
metric soft mode when the vortex is situated at the cen- 
ter of the system. Switching can occur only as a result of 
nonlinear mixing between the radially symmetric mode 
and non-symmetric vortex modes which do interact with 
the spatially uniform alternating external field. This may 
take place in large systems where the motion of the vortex 
is frozen Jig]. However, in a relatively small, finite system 
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with free boundary conditions the vortex is attracted by 
its image and moves along an unwinding spiral trajectory 
(see e.g. $7\ ) towards the boundary. Our numerical ex- 
periments do show that switching events in general occur 
only when vortex is at a finite distance from the center. 
The vortex center motion is very slow (with a frequency 
~ 1/L 2 ). So we can consider the switching process with 
fixed vortex position, say X — Rcosx,Y = Rsm\- In- 
serting Eqs (|J) and (^) into Eq. (||) and assuming that 
the vortex is far from the boundaries (R <C L), we find 
that the effective interaction of the in-plane ac magnetic 
field with the soft mode: 



\ 7 t) = h [ cii ipi sin(cji) — h^(a 2 + bmf) cos(ujt) 



where a ; = J2n 
'■v' 



x 2 



(8) 



(1C Hii y, (l = 1,2), and b = 



2L 2 



- (Cfii) ■ An effective interaction of the 



same form can be obtained by taking into account the 
fact that the vortex structure is velocity dependent |lC| ] 
and symmetric about the direction of the vortex motion. 
The constant phase shift (q\) plays no essential role and 
was omitted. 

From Eqs (0) and (||) we find that in the soft mode 
approach the dynamics is governed by 



mi 



h (ai sinwi — 02^1 coswi) 
-7(/ziTOi + Ami), 
ipl = Himx + Amf — T^i — hbm\ cosut. 



(9) 



Here the nonlinear terms m™, (n > 3) and m\ip\ have 
been neglected. An example of the core dynamics based 
on Eqs (||) is presented in Fig. 3. These results are in 
good agreement with those from the numerical integra- 
tion of the full Landau-Lifshitz equations (Q) . 

To clarify the physical meaning it is convenient to 
write the set of equations (|9|) as a single equation for 
m\. Near the threshold ( \ni\ -C 1), in the limit of 
small damping (7 -C 1), and for a not too strong 
amplitude of the external magnetic field ( a 2 h < 1) 
we obtain from Eqs (^J) an effective equation for mi 
rrii + J nix + /ii mi + Ami + h(ai li — bm±) cosu>t = 0. 
Thus the vortex core dynamics is analogous to the dy- 
namics of a particle in a double-well potential under the 
action of direct and parametric forces. In order to esti- 
mate h cr we use a heuristic approach proposed by Moon 
fL8| . Here the switching occurs when the particle reaches 
the maximum velocity on the homoclinic orbit. Consid- 
ering the particle motion in the potential well centered at 
mi = p ^JJli and using multiple-scale analysis for small 
7, h and \uj 2 — uJq\ <§; 1, where uj — ^/2|/ii| is the fre- 
quency of harmonic oscillations near the bottom of the 
well, we find m\ ~ yj/Ii | (p + M{u>) cos(wt)), with 



M 



0^0 M' 



1 2 2 
+ 7 U) 



h 2 

2A— (ai lu — b\ 



(10) 
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The maximum velocity on the homoclinic orbit is 
J|^-.The switching occurs when wM(u>) — a^, where 
a sa 1 is an empirical parameter. From Eq. ( J10| ) we find 
h cr (u>) in the form 



f2A |Q (ai-Jttl - bp) 




n 2 



3 a 2 
8^2 



^2 



(ii) 



where Q = lu/loq. The condition (|Tl]) is in agreement with 
the results of numerical simulations: the function h cr {u>) 
has minima near the frequency of the soft mode u> ~ ±o->o- 
It is highly asymmetric ( Fig. 4). For small amplitudes 
of the counter-clockwise rotating magnetic field (uj > 0), 
the switching condition h > h cr can be fulfilled only 
for the vortex with initial polarization p = — 1, while for 
a clockwise rotating field the condition can be fulfilled 
for the vortex with opposite polarization p = 1. Quali- 
tatively the same results may be obtained by using the 
Melnikov function approach [|l9| . 

In conclusion, we have shown that the polarization of 
out-of- plane vortices in easy-plane ferromagnets can be 
changed by applying an ac magnetic field. Flips occur 
more easily when the polarization of the vortex is anti- 
parallel to the angular velocity i2. It can take place only 
in discrete systems, and they are uni-directional events. 
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Figure 2: Threshold value of the held amplitude h cr 
versus u> obtained from the numerical integration of the 
full Landau-Lifshitz equations. 
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Figure 1: Switching from the state with positive po- 
larization to the state with negative polarization of the 
out-of-plane vortex due to a clockwise rotating magnetic 
field with the frequency uj = —0.1. The damping con- 
stant 7 = 0.002. The lower curve shows the time evolu- 
tion of the magnetization of the inner shell when the field 
amplitude h = 3 x 10 -3 is above the threshold value, the 
straight line corresponds to h = 10~ 3 (the amplitude of 
oscillations in this case is very small w 0.005 and there- 
fore they are not seen in the figure) 



Figure 3: Time evolution of the vortex core magneti- 
zation in the presence of the clockwise ( upper curve) and 
counter-clockwise (lower curve) rotating magnetic field 
based on Eqs (^|). The parameters are A = ai = a-2 — 
b=l,[i 1 = -0.1, uj = ±0.1, h = 0.04,7 = 0.01. Initial 
polarization p = +1. 
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Figure 4: Threshold value of the field amplitude h cr 
versus ui, as given by Eq. (pd|). The parameters used are 
01 = b,Ab'/a = 0.01, w = 0.08, 7 = 0.002. 
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